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Abstract 

To offer an insight into the rapidly developing theory of fractional 
diffusion processes we describe in some detail three topics of current 
interest: (i) the well-scaled passage to the limit from continuous time 
random walk under power law assumptions to space-time fractional 
diffusion, (ii) the asymptotic universality of the Mittag-Leffler waiting 
time law in time-fractional processes, (iii) our method of parametric 
subordination for generating particle trajectories. 
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1 Introduction 



The field of fractional (more generally anomalous) diffusion processes in 
recent decades has won more and more interest in applications in the 
sciences, in physics and chemistry, and even in finance. Here we will give 
some insight into this rapidly developing field. Instead of trying to cover the 
whole spectrum of this field wc will focus on three contrasting aspects to 
which we ourselves have affinity by our research, thereby trying to keep this 
paper as self-contained as possible. We hope not only to whet the appetite 
among people not fully familiar with the subject but also to propagate our 
methodological viewpoints in the fractional diffusion community. 

Viewing fractional diffusion processes as generalization of the familiar 
classical diffusion process, we begin by considering the Cauchy problem for 
the classical diffusion equation 

du(x,t) d^u(x,t) , „, , „ „ 

— ^ = ' u{x,0+) = f{x), xe-R,t>0. (1.1) 

In fractional diffusion equations (see Section 2 for explanations) the 
differentiations with respect to t and x are replaced by differentiations of 
non-integer order. For equation (1.1) it is well known that its solution u{x, t) 
has the essential properties we expect from a diffusion process, that is a 
process of re-distribution in space x and time t. Considering u{x, t) as the 
spatial density of an extensive quantity, e.g. mass, charge, or probability, 
we have (a), (b) and (c): 

(a) conservation of the total quantity: J^^u{x,t) dx = J^^u{x,0'^) dx, 
Vt > 0. 

(b) preservation of non-negativity: u{x,Q~^) > 0, Vx G R implies u{x,t) > 
0, Vx G R, t > 0. 

(c) Another essential characteristic of problem (1.1) concerns the law of 
spreading (or dispersion) of the quantity. With the special initial condition 
u{x, 0"*") = S{x) (the Dirac generalized function), the variance grows linearly 
in time, that is cJ^(t) := x"^ u{x, t) dx = 2t. More generally, in a classical 
diffusion process the variance, which is a natural and common quadratic 
measure of the spread of a diffusing substance, grows linearly in time, that 
is, if we allow a drift, we have o"^(t) := x'^ t) — m{t)] dx ^ Ct as 
t ^ oo with m{t) := Jl^^ xu{x,t) dx, for some constant C > 0. 

The above properties (a), (b) and (c) are indeed shared by many processes 
governed by second-order linear parabolic equations. Usurping the term 
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diffusion for processes having properties (a) and (b) but not necessarily 
(c), we follow the custom to call processes of anomalous diffusion those in 
which, for initial condition u{x,0~^) = 5{x), the variance does not exhibit 
essentially linear grow with i — > oo. Among these processes we single out 
the sub-diffusive ones for which the variance grows (for large t) more slowly 
than linearly, and the super- diffusive for which it grows (for large t) faster 
than linearly, or even does not exist (i.e is infinite). Focusing our attention 
to the space-time fractional diffusion equation (i.e. a generalization of the 
classical diffusion equation (1.1) via suitable pseudo-differential operators 
interpreted as time and space derivatives of fractional order), we will discuss 
the construction and properties of its fundamental solution (obtained for 
f{x) = 5{x)) and its approximation by continuous time random walk. This 
generalized diffusion equation has found wide interest among researchers 
in recent 20 years. After a general survey of basic facts we will go into 
details of three distinct but related topics. As this paper cannot be a 
substitute for an extensive monograph, our presentation will naturally be 
biased by our and our close co-workers' contributions. We will meet the two 
complementary aspects of diffusive processes. The first is the macroscopic 
aspect: the structure of the fundamental solutions in dependence on space 
X and time t, in particular their scaling properties and asymptotics. Here 
u{x, t) is viewed as the density with respect to x evolving in t. The second is 
the microscopic aspect: the structure of the trajectories (paths) of particles 
subject to such process. Here u{x, t) is viewed as the sojourn probability 
density with respect to x evolving in t. 

The rough structure of our paper is as follows. In Section 2 we will give 
a survey of the space-time fractional diffusion equation and the essential 
properties of its fundamental solution. Sections 3 and 4 are devoted to topic 
(i): continuous time random walk (CTRW), Section 5 to topic (ii): power 
laws and well- scaled passage to the diffusion limit, Section 6 to topic (iii): 
our method of parametric subordination for exact simulation of trajectories. 
Finally, in Section 7, we will draw some conclusions. 

Of course, there are more significant advances than we can report here. Let 
us only mention processes with distributed orders of fractional differentiation 
and multi-dimensional processes and cite, e.g. the papers [H El El SS] and 
the fundamental monograph by Meerschaert & Scheffler |33]. Much work 
has been done in numerical analysis of difference schemes for calculating 
densities of fractional diffusion processes, see e.g. [22] . We apologize to all 
authors who feel that we appreciate their work insufficiently. Throughout 
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this paper we will make liberal use of generalized functions in the sense of 
Gel'fand and Shilov [10], so avoiding the cumbersome notations of measures 
and functionals, and usually we will assume the occurring (generalized) 
functions so well-behaved that our manipulations are allowed. By this we 
hope our presentation to be welcome to researchers in applications as well as 
inspiring for pure mathematicians looking around outside the ivory tower. 



2 The space-time fractional diffusion 

We begin by considering the Cauchy problem for the (spatially one- 
dimensional) space-time fractional diffusion equation 

tD^ u{x,t) = ^D^u{x,t) , u{x,0)=6{x), xeK, t>0, (2.1) 

where a, 9, (3 are real parameters restricted to the ranges 

0<a<2, |6l| <min{a,2-a}, < /? < 1 . (2.2) 

Here t-D=f denotes the Caputo-Dzherbashyan fractional derivative of order 
P, acting on a sufficiently well-behaved function f{t) of the time variable t, 

and xDq denotes the Riesz-Feller fractional derivative of order a and 
skewness 9, acting on a sufficiently well-behaved function g{x) of the space 
variable x, 



.D^gix) = + jsin [{a + 9)^] 



g{x + x') - g{x)_ 



+ sin 







poo 


(a - ( 










Jo 



(2.4a) 

g{x - x') - g{x) ^^,| ^ Q ^ ^ ^ 2 , \9\ < min{a, 2-a} . 



In the symmetric case = 0, which later will be of our main interest, formula 
(2.4a) simplifies to 

Vil + a) . (a-K\ f°° g(x + x')-2g(x)+g(x-x') 

sm I / : 



xDn f(x) = sin ( — I / ^- dx' . 

(2.4'a) 

The above representations of the space fractional derivatives are based on a 
suitable regularization of hyper-singular integrals. 
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In the limits /3 = 1 and a = 2 (so 9 = 0) we recover the first time derivative 

df{t) , , , , . . d^gix) . , 

— ; — and the second space derivative — t"^, respectively. 
dt ax^ 

These representations mirror the fact that time-fractional (for < /? < 
1) processes are processes with long memory (see for this also Section 5) 
whereas space fractional (for < a < 2) are processes with spatial long- 
range interactions. For more information on these operators we refer to 
[ia[l5l|26l[28l|38llin]. 

Let us note that the solution u{x,t) of the Cauchy problem (2.1), known as 
its Green function or fundamental solution, can be viewed as the density of 
an extensive quantity or as the probability density in the spatial variable x, 
evolving in time t. In the case a = 2 (hence ^ = 0) and /? = 1 we recover 
the standard diffusion equation for which the fundamental solution is the 
Gaussian density with variance cr^ = 2t. 

Writing, with Re[s] > ctq, k G R, the transforms of Laplace and Fourier as 

C{fity,s} = f{s):= / e-'^f{t)dt, 
Jo 

/ + 00 , 
e"^^g{x)dx, 
-oo 

the corresponding transforms of tD^ f{t) and xDQg{x) are 

>C { tD^ f{t)] = sl^J{s)-s(^-^ /(O) , (2.36) 
J^{^D'^g{x)] = -|Ac|«i^signK^(^) ^ _|^| a (sign k) ^ 7r/2 ^ (2.46) 



We will freely use the convolution theorems pertinent to these transforms, 
defining for generic functions: (/i * f2){t) = /jq t] ~ /2(^') ^^'j * > 0, 
(51 * 92){x) = /(_oo,+cx)) - x') g2{x') dx' , x G R, and the convolution 
powers f*^{t) and g*^{x) as n- fold convolutions ( n > 0). Note that f*^{t) = 
6{t), g*^{x) = 6{x). For mathematical details we cite [HI [38] on the Caputo- 
Dzherbashyan derivative and [40j on the Feller potentials. For the general 
theory of pseudo-differential operators and related Markov processes the 
interested reader is referred to the excellent volumes by Jacob [23] . 

Let us here recall the representation in the Fourier-Laplace domain of the 
(fundamental) solution of (2.1). Using 6{k) = 1 we have from (2.1) 

s ^ S(k, s) - s /5 - 1 = - 1^1 « i ^ sign K ^) ^ (2.5) 
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hence explicitly 



u{k, s) 



(2.6) 



s/3 + |^|a^(9signK ' 



For explicit expressions and plots of the fundamental solution of (2.1) in the 
space-time domain we refer the reader to 128] . There, starting from the fact 
that the Fourier transform u{k, t) can be written as a Mittag-LefSer function 
with complex argument, the authors have derived a Mellin-Barnes integral 
representation of n(x, t) with which they have proved the non-negativity of 
the solution for values of the parameters {a, 0, /3} in the range (2.2) and 
analyzed the evolution in time of its moments. In particular for {0 < a < 
2, /3 = 1} we obtain the densities of the stable processes of order a and 
skewness 0. The representation of n(x, t) in terms of Fox //-functions can 
be found in [29]. We note, however, that the solution of the space-time 
fractional diffusion equation (2.1) and its variants has been investigated by 
several authors as pointed out in the bibliography in [28j: here we refer to 
some of them, [21 [3l [311 [32l [3l], where the connection with the CTRW was 
also pointed out. 

Henceforth our attention, if not said explicitly otherwise, will be focussed 
on the symmetric case = 0. In this case u{x,t) is an even function of x 
and we get from (2.6) the Fourier-Laplace representation 



which allows us demonstration of the conservation property (a), namely 
f^^u{x,t) dx = 1 for all t > 0, and calculation of the variance cT^(t) = 
(x^(t)) (the quadratic measure of the spread as function of t) of a diffusing 
particle. 

From (2.7), more generally already from (2.6), we find by aid of well-known 
properties of the Fourier transform u{0,s) = 1/s, hence u{x,t) dx = 
u{0, t) = I and so /J*"^ u{x, t) dx = 1. 
For the variance we find 



u{k, s) 



(2.7) 




a 




(2.8) 



Using the Mittag-Leffier function 




(2.9) 
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see e.g. |141 138 j . we find by Laplace inversion the convergent series 



r(i + /3) r(i + 2/3) 

from which for t > we get 

( 2t , a = 2 , (3 = 1 (normal diffusion) , 

12+/3 
^^^^p^ ' a = 2 , < /3 < 1 (sub-diffusion) , (2.10) 
oo , 0<Q!<2,0</3<1 (super-diffusion) . 

Admitting that these calculations are formal we leave the task of strict 
justification to pure mathematicians. 

3 The continuous time random walk 

In the Sixties and Seventies of the past century Montroll and Weiss and Scher 
(just to cite these pioneers) published a series of papers for modelling rather 
general diffusion processes by random walks on lattices, see e.g. |35^ l36j. and 
the book by Weiss [48j with references therein. Initiated by their activity 
the concept of continuous time random walk became popular in physics. 
CTRWs are good phenomenological models for several types of diffusion, in 
particular the microscopic aspects of a particle jumping from point to point 
admitting waiting times between jumps. Allowing all space instead of only a 
lattice, a CTRW can be viewed as a compound renewal process or a renewal 
process with reward (see [7J), or a random walk subordinated to a renewal 
process. Let us recall the basic notions of the CTRW theory. 

A CTRW is generated by a sequence of independent identically distributed 
(iid) positive random waiting times Ti , T2 , T3 , . . . , each having the same 
probability density function <j){t) , t > , and a sequence of iid random 
jumps Xi, X2, X3, . . . , in R , each having the same probability density w{x) , 
x € R . Setting to = , t„ = Ti -|- T2 + . . . T„ for n S N , the wandering 
particle makes a jump of length in instant t„, so that its position is 
xo = for < i < Ti = ti , and Xn = Xi + X2 + . . . Xn , for t„ < i < tn+i ■ 
We require the distribution of the waiting times and that of the jumps to 
be independent of each other. We allow the probability densities and w 
to be generalized functions in the sense of Gel'fand and Shilov [TOJ. 

Natural probabilistic arguments lead us to the integral equation for the 
probability density p{x, t) (a density with respect to the variable x) of the 
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particle being in point x at instant t , see e.g. [191 EOl Hll SSI E], 



w{x — x') p{x' , t') dx' 



dt' , (3.1) 



in which the survival function ^{t) = (j){t')dt' denotes the probabihty 
that at instant t the particle is still sitting in its starting position x = . 
Clearly, (3.1) satisfies the initial condition p{x,0) = 5{x), and p in place of 
u has the properties (a) and (b) of Section 1. In the Fourier-Laplace domain 
Eq. (3.1) reads as 

p{k, s) = $(s) + w{k) (j){s)p{K, s) , (3.2) 
and using ^'(s) = (1 — (j){s))/s , explicitly 

= (3.3) 

This representation is known as the the Montroll- Weiss equation., so named 
after the authors of |36j, who derived it in 1965 as the basic equation for 
the CTRW. By inverting the transforms one can find the evolution of the 
sojourn density p{x,t) for time t running from zero to oo. In fact, recalling 
that |u;(k)| < 1 and \(p{s)\ < 1, if k 7^ and s / 0, Eq. (3.3) becomes 

00 00 

f{K,s) = $(S) Y,[^{s)w{K)r = J2^n{s)Wn{K), (3.4) 
n=0 n=0 

and we promptly obtain the series representation of the continuous time 
random walk, see e.g. [7] (Ch. 8, Eq. (4)) or ^ (Eq.(2.101)), 

oo oo 
p{x, t) = Y^ Vn{t) Wn{x) = ^{t) 8{x) + ^ V„(t) Wn{x) , (3.5) 
71=0 n=l 

where the functions Vn{t) and Wn{x) are obtained by iterated convolutions 
in time and in space, Vn{t) = (^' * </)*"')(t), and Wn{x) = {w*"'){x), 
respectively. In particular, vo{t) = * 5){t) = ^{t), vi{t) = * (p){t), 
wo{x) = 5{x), wi{x) = w{x). In the R.H.S of Eq (3.5) we have isolated 
the first singular term related to the initial condition p{x,0) = ^(0) 6{x) = 
5{x). The representation (3.5) can be found without detour over (3.1) by 
direct probabilistic reasoning and transparently exhibits the CTRW as a 
subordination of a random walk to a renewal process: it can be used as 
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starting point to derive the Montroll- Weiss equation, as it was originally 
recognized by Montroll and Weiss [36]. 

With the special choice (j){t) = me""**, m > 0, equation (3.1) describes the 
compound Poisson process. It reduces after some manipulations (best carried 
out in the Fourier-Laplace domain) to the Kolmogorov- Feller equation: 

— p{x,t) = —mp{x,t) + m / w{x — x') p{x\t) dx' . (3.6) 

J —oo 

Then from (3.5) we obtain the series representation 

°° /' 4.\k 

PM = Y,^e-"^'wu{x). (3.7) 

A:=0 

Note that only in this case the corresponding stochastic process is 
Markovian. 



4 Relevance of power laws and well-scaled passage 
to the diffusion limit 

In this section we work out the effect the power laws for the distributions of 
jumps and waiting times have on the limiting properties of continuous time 
random walks. By power law we usually mean a law of decaying at infinity 
like a power (with negative exponent) of the relevant independent variable. 
There is a vast amount of literature on such laws. Let us recommend just 
two items, namely ^37] and [IS]. What we are going to show now is that 
appropriate power laws for jumps and waiting times are microscopic models 
for fractional diffusion processes. 

To be sufficiently general we introduce the cumulative functions for waiting 
times and jumps. With our densities (j){t) and w{x) we define 



t X 



-OO < X < oo. (4.1) 



These functions are non-decreasing, <I>(0) = 1^(— oo) = 0, $(00) = VF(oo) = 1. 
As they may have points of discontinuity we agree on the provision that 
equations in which they occur are meant to hold at points of continuity. 
The notion of power law concerns the behaviour of <I>(t) and W{x) for large 
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t and large |a;|, roughly the mode of decrease of the tails 1 — <I>(t) = ^{t) 
near t = oo and W{x) near x = — oo, 1 — W{x) near x = +00, like a power 
of t or |x| with a negative exponent. 

Not wanting to overload our presentation we assume spatial symmetry (6 = 
0) with respect to x = and avoid decoration of asymptotic behaviours 
with slowly varying functions. For such neglected generalizations we refer 
to [m |T7]. The parameters a and /3 of equation (2.1) play essential roles 
yielding power laws in the strict sense ifO<a<2,0</3<l, but still 
formal analogies in the extreme cases a = 2, f3 = 1. 

The question of interest is the long-time, wide-space behaviour of a CTRW 
under power law assumptions for waiting times and jumps, i.e. the 
appearance of such CTRW when observed after long time and from far 
away. To bring the distant future and the far-away space into near sight 
we multiply time intervals by a small positive factor r, space intervals by a 
small positive factor h, so making large intervals numerically of moderate 
size, intervals of moderate size numerically small. Essentially this means 
changing the units of time and space from 1 to 1/r and 1/h, respectively. 
We then obtain the asymptotic behaviour by sending r and h to zero, in a 
specially combined way, namely under the requirement of a scaling relation, 
honouring which we call the whole procedure "well-scaled passage to the 
diffusion limif. In this limit we will obtain a process obeying the space- 
time fractional diffusion equation (2.1) with 6 = , a fact that justifies 
the CTRW approach. Conversely, we can consider a CTRW as a model "in 
the small" of a space-time fractional diffusion process. For generalization to 
skewed processes {9 7^ 0) see [17] . 

As we carry out the essential work in the Fourier-Laplace domain we use 
the fact that the behaviour of functions f{t), g{x) in the infinite is mirrored 
in that of their transforms /(s), ^(k) for s ^ 0"*", k ^ 0. The lemmata 
we need are provided by the Tauber-Karamata theory for which we refer 
to [1] and [S], they can also be distilled from the Gnedenko theorem on the 
domains of attraction of stable probability laws, see [TI]. See [13] and [T7] 
for more general versions. What we need is contained in the following two 
lemmata which are simplified versions of more general facts. For the reader's 
convenience we give proofs in the Appendix. 
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MASTER LEMMA 1: 

AssumeW{x) increasing: W{—oo)=0,W{+oo) = l, 

symmetry: J dW{x') = J dW{x') for x > 0, and either (a) or (b): 

{—oo,—x) (a;,+oo) 

(a) := J x'^dW{x) < oo, labelled as a = 2 , 

(— oo,+oo) 

(b) J dW{x') ~ 6a-^x-" for x +oo, a G (0,2) and 6 > 0. 

Then we have the asymptotics I — w{k,) ~ fJ'\K\°' for k ^ with = in 
case (a) and /i = 67r/[r(a + 1) sin(a7r/2)] in case (b). 
MASTER LEMMA 2: 

Assume increasing: $(0) = 0, $(+cxd) = 1, and either (A) or (B): 

(A) p := f td^{t) < oo, labelled as (3 = 1, 

(0,+oo) 

(B) m{t)= J d^it)'^c(3-H-'^ fort -^oo, (3 £ (0,1) and oO. 

(t,oo) ^ 

Then we have the asymptotics 1 — (p^s) ~ Xs^ for < s — > with X = p in 
case (A) and X = C7r/[r(/3+ 1) sin(/37r)] in case (B). 

Assuming the conditions of these two lemmata satisfied we are ready for 
passing to the diffusion limit. We multiply the jumps by a factor ^ > 0, 
the waiting times by a factor r > 0. So we get a transformed random 

n n 

walk Sn{h,T) = hX}^ with jump instants tnir) = ^ tT^ that we now 

k=l k=l 
investigate with the aim of passing to the limit h ^ 0, t ^ under a 

scaling relation between h and r yet to be established. As it is convenient 

to work in the Fourier-Laplace domain we note that the density (prit) of 

the reduced waiting times rT/j and the density Wh{x) of the reduced jumps 

hXk are ^r(i) = (t>{i/T)/T, t > ; Whix) = w{x/h)/h, — oo < x < oo. The 

corresponding transforms are simply (pris) = 4>{st), Whin) = w{Kh). We are 

interested in the sojourn probability density Ph,Tix, t) of the particle subject 

to the transformed random walk. In analogy to the Montroll- Weiss equation 

(3.3) we get 

|,,(„,.) = 1^_J__ = '-^>'' . (4.2) 

S 1 — Wh{K)4>r{s) * 1 — w{hK)(j){sT) 

Considering now s and k fixed and 7^ we find for ^ — 0, r ^ from the 
Master Lemmata (replacing there k by nh, s by sr ) by a trivial calculation, 
omitting asymptotically negligible terms, the asymptotics (4.3) with (4.4). 

PhA^^ 5) ~ + xirsf ~ r{h, t)|/«|" + sP ' ^ ' 
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So we see that, for every fixed real k 7^ and positive s, 



as /i and r tend to zero under the scaling relation r{h,T) = 1. Comparing 
with (2.7) we recognize here u{n, s) as the combined Fourier-Laplace 
transform of the solution to the Cauchy problem (2.1) with 9 = 0. 
Invoking now the continuity theorems of probability theory (compare [8]), 
bypassing some analytical subtleties, we see that the time-parametrized 
sojourn probability density converges weakly (or in law) to the solution of 
the Cauchy problem (2.1) with 9 = 0. This weak convergence can be taken 
as justification for approximate simulation of trajectories (particle paths) by 
CTRW's with power law jumps and waiting times so chosen that routines 
are available for producing the needed random numbers (see e.g. [13j). 

5 The Mittag-LefHer waiting time law and time- 
fractional processes 

We now offer another view to the transition from CTRW (under power law 
assumptions) to fractional diffusion, separating the temporal and spatial 
limiting procedures, thereby getting among other matters the time- fractional 
CTRW discussed in the pioneering paper [21] . 

Turning our attention to time-fractional processes we present in a condensed 
way some results from our papers [181 US ED]- We will see the importance 
of the Mittag-Leffler waiting time density (j)^^{t) and the corresponding 
survival function '^^^^{t) with their Laplace transforms displayed here: 

0*^^(0 = -{d/dt)Ep{-tP) , ^"'^{s) = , (5.1) 

^^\t) = Ep{-tP), $^«(.) = ^f^, (5.2) 

with Ep defined in (2.9). Throughout we assume < (3 < 1. We recall that 
for /? = 1 we recover (j)'^^{t) = exp(-t), ^^^^(t) = exp(-t). 

The importance of these functions cannot be overestimated, they also play 
an essential role in the theory of fractional relaxation (see e.g. |14j). The 
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relevance of the Mittag-Leffler waiting time law for time-fractional CTRW 
has been put in bright light by Hilfer and Anton in 1995 [21] • Fulger, Scalas 
and Germano [9] pay special attention to its use as waiting time law in 
CTRW simulation . In the Sixties of the past century it has been found by 
Gnedenko and Kovalenko [12] as a limiting law in thinning (rarefaction) of 
a renewal process under the power law assumptions of our Master Lemma 
2, but they only gave its Laplace transform without identifying it as a 
function of Mittag-Leffler type. We will show, without invoking the concept 
of thinning, that it is a universal limiting law for long-time behaviour of a 
renewal process under a power law regime. 

For a general CTRW we will show how via the concept of the memory 
function we can separate the passages of the scaling factors r and h (of 
the preceding Section 4) to zero, thus avoiding the simultaneous use of the 
continuity theorems for the transforms of Laplace and Fourier. 



5.1 Manipulations: rescaling and respeeding 

To introduce the memory function H{t) and explain its meaning we recall 
the CTRW theory of Section 3, in particular equations (3.1) and (3.3). We 
need this general theory already for embedding into it the renewal theory. 
In fact: we can view a pure renewal process as a special type of CTRW, 
namely one in which all jumps are of fixed size 1 and the position x of 
the wandering particle in space plays the role of the counting number n of 
the renewal process. However, we prefer to continue working in the general 
CTRW context. Introducing formally in the Laplace domain the auxiliary 
function 

= izilfl = Sfl , hence 0(s) = (5.3) 

scl){s) (Pis) 1 + sHis) 

and assuming that its Laplace inverse H{t) exists, we get, following [30], in 
the Fourier-Laplace domain the equation 

H(s) sp{k, s) - I = [w{k) - 1] p{k, s) , (5.4) 

and in the space-time domain the generalized Kolmogorov-Feller equation 
/■* d /■+°° 

/ H(t -t') —p(x,t')dt' = -p(x,t) + w(x - x')p(x',t)dx', (5.5) 

Jo of 
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with p{x,0) = S{x). With the special choice of the power law memory 
function 

i^^^(t)= rxl^' o</^<i' (5.6) 

{5{t), if (3 = 1, 

whose Laplace transform is 

H^^{s) = sl^-^ , 0</3<l, (5.7) 

wc have the Mittag-Leffler waiting time law given by Eqs (5.1) and (5.2). 
In the extremal case /? = 1 this reduces to the exponential waiting time 
law (j){t) = exp(— i), '^{t) = exp(— i), and we obtain in the Fourier- Laplace 
domain 

sp{k, s) — 1 = [w{k) — 1] p{k, s) , (5-8) 
in the space-time domain the classical Kolmogorov-Feller equation 

—p{x,t) = —p{x,t)+ I w{x — x')p{x',t)dx' , p{x,0) = S{x). (5.9) 

For < /3 < 1 we have the time-fractional Kolmogorov-Feller equation 

/+0O 
w{x — x') p{x' , t) dx' , p{x, 0"*") = 5{x) . (5.10) 
-oo 

Let us now consider two types of manipulations on the CTRW: 
A: rescaling, B: respeeding. 

(A) means, as in Section 4, change of the unit of time. With the positive 
scaling factor T we replace the waiting time T hy tT , again intending r << 1. 
In a moderate span of time we will so have a large number of jump events. 
Again we get the rescaled waiting time density and its corresponding Laplace 
transform as (j)Tit) = (p{t/T)/T, (pris) = <P{ts)- By decorating also the 
density p with an index r we obtain this rescaled CTRW integral equation 
in the Fourier-Laplace domain as 

Hr{s) sp^{k, s) - 1 = [w{k) - 1] p^{k, s) , (5.11) 

where _ 

H,{s) = ^ ~ f , hence ^(rs) = 1--. (5.12) 

s(p{TS) 1-\-sHt{s) 
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Remark: Note that in this Section 5 the position of the indices at the 
density p and their meaning are convenient but different from those in the 
other Sections. 

d 

(B) means multiplying the quantity representing —p(x,t) by a factor 

1/a, where a > is the respeeding factor, a > I means acceleration, 
< a < 1 means deceleration. In the Fourier-Laplace representation this 
means multiplying the RHS of Eq. (5.11) by the factor a since the expression 
sp{k,s) — 1 , in view of p{x,0) = 5{x) and 5{k) = 1, corresponds to 

We now consider the procedures (A) and (B) in their combination so that 
in the transformed domain the rescaled and respeeded process has the form 



(5.13) 



Clearly, the two manipulations can be discussed separately: the choice {r > 
0, a = 1} means pure rescaling, the choice {r = I, a > 0} means pure 
respeeding of the original process. In the special case r = 1 we only respeed 
the system; if < r <C 1 we can counteract the compression effected by 
rescaling to again obtain a moderate number of events in a moderate span 
of time by respeeding (decelerating) with < a <C 1. These vague notions 
will become clear as soon as we consider power law waiting times. Defining 
now _ _ 

Hr{s) 1 - <PiTs) 



as(p{Ts) 

we finally get, in analogy to (5.11), the equation 



SPr 



1 



= [w{k) - 1] p {k, s) . 



(5.14) 



(5.15) 



What is the combined effect of rescaling and respeeding on the waiting time 
density? In analogy to (5.3) and taking account of (5.14) we find 

1 1 



</'r,a(s) 



1 + sHrJs) 



1 + S 



1 - (i){Ts) 



(5.16) 



as (t>{Ts) 

and so, for the deformation of the waiting time density, the essential formula 



TS 



1 - (1 - a)(^(Ts) 



(5.17) 
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5.2 Asymptotic universality of the Mittag-LefHer waiting 
time law under power law regime 

We now recall the MASTER LEMMA 2 of Section 4 and assume the 
conditions stipulated there. By using the statements of this lemma, taking 

a = XT^, (5.18) 

fixing s as required by the continuity theorem of probability for Laplace 
transforms, the asymptotics (f){s) ~ 1 — As^ for < r ^ implies 

Ar^ [1 - XtI^sI^ + o{tI^sI^)] 1 
'^^'^^'^ " 1 - (1 - At/?) [1 - Xt^sP + o{rPs(^)] ^ 1+^ ' ^^'^^^ 

corresponding to the density 4>^^{t). This formula expresses the 
asymptotic universality of the Mittag-LefHer waiting time law that 

includes the exponential law for (3 = 1. It says that our general power law 
waiting time density is gradually deformed into the Mittag-Leffler waiting 
time density. 

Remark: Let us stress here the distinguished character of the Mittag-LefHer 
waiting time density (f)^^{t) defined in (5.1). It is easy to prove the identity 

4>f^^{s) = 4>^^{Ts/a^"^) for all r > 0, a>0, (5.20) 

that states the self-similarity of the combined operation rescaling-respeeding 
for the Mittag-Leffler waiting time density. In fact, (5.20) implies (p^^ait) = 
(l)^^(t/c)/c with c = t/o^/^ , which means replacing the random waiting 
time by cT^^. As a consequence, choosing a = t^, we have 

0^,^,(s) = 0^^(s) for all r>0. (5.21) 

Hence the Mittag-Leffler waiting time density is invariant against combined 
rescaling with r and respeeding with a = t^. Observing (5.19) we can say 
that 4>'^^{t) is a r — > attractor for any power law waiting time (compare 
Master Lemma 2) with 

^{t)^^t~'^, 0</3<l, OO, (5.22) 

y 

under combined rescaling with r and respeeding with a = Xt^. This 
attraction property of the Mittag-Leffler waiting time distribution with 
respect to power law waiting times (with < (3 < l)isa kind of 
analogy to the attraction of sums of power law jump distributions by stable 
distributions. 
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5.3 Diffusion limit in space 

We can obtain from (5.5) the fractional Kolmogorov- Feller equation (5.10) 
for time fractional CTRW by direct insertion of the Mittag-Leffler memory 
function into the equation or, as the previous considerations show, by 
manipulating it under power law assumption for the waiting time and 
passing to the limit. We have not yet operated on the jumps. To do this now, 
we assume the conditions of Master Lemma 1 to hold. Then, by another 
respeeding, in fact an acceleration (that we earlier had carried out in jl9j). 
we will arrive at diffusion processes fractional in space. We have three 
choices: 

(A) : dijfusion limit in space only, for general waiting time, 

(B) : diffusion limit in space only, for Mittag-Leffler waiting time, 

(C) : joint limit in time and space (with power laws in both) with scaling 
relation. 

Note that (B) is just a special case of (A) but of particular relevance. In 
all three cases we rescale the jump density by a factor h > 0, replacing the 
random jumps X by hX. This means changing the unit of measurement in 
space from 1 to l//i, with < /i ^ 1, so bringing into near sight the far-away 
space. The rescaled jump density is Wh{x) = w{x/h)/h, corresponding to 

Whin) = WfJlHi). 

Choice (A): Starting from the Eq. (5.4), the Fourier-Laplace representa- 
tion of the CTRW equation, without special assumption on the waiting time 
density, we accelerate the spatially rescaled process by the respeeding factor 
l/{fih"), arriving at the equation (using qh as new dependent variable) 



His) 



^—L—q^{K,s). (5.23) 



Then, fixing n as required by the continuity theorem of probability theory 
for Fourier transforms, and sending h to zero we get, noting that by Master 
Lemma 1 [w{hK) — l]/(/i/i") — > — and writing u in place of qq, 

H{s) su{k,s)-1 = -|k|"u(k,s) , (5.24) 

still with H{s) as in (5.3). Here — is the symbol of the Riesz pseudo- 
differential operator xDq (known as the Riesz fractional derivative of order 
a) obtained from the Riesz- Feller fractional derivative for 9 = Q, see (2.1) 
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and (2.4). We thus arrive at the integro-pseudo-differential equation 

/ H{t - t') -^uix, t') dt' = ^D^ u{x, t), 0<a<2, u{x, 0) = d{x) . 
Jo 

(5.25) 

Comments: By this rescahng and acceleration the jumps become smaller 
and smaller, their number in a given span of time larger and larger, the 
waiting times between jumps smaller and smaller. In the limit there are 
no waiting times anymore, the original waiting time density (f){t) is now 
only spiritual, but still determines via H{t) the memory of the process. Eq. 
(5.25) offers a great variety of diffusion processes with memory depending 
on the choice of the function H{t). 

Choice (B): Inserting in (5.25) the Mittag-Leffler memory function (5.6), 
we immediately get the space-time fractional diffusion equation (2.1) with 

e = o. 

Choice (C): Assuming the conditions of both Master lemmata fulfilled, 
rescaling as described the waiting times and the jumps by factors r and h, 
starting from (5.13), decelerating by a factor Ar^ in time, then accelerating 
for space by a factor l/{ij,h°'), we obtain, by fixing s and k, the equation 

Hris) [sPr,a{T,h)i'^, S) " ij = a{T, h) [wh{K) - 1] Pr,a{T,h)i^i^ s) , 

with Wh{K) = w{hK), a{T, h) = Xt^^ / {iJ,h°') and 

H^^s)=^^^3^^Xt^s^-\ for T^O. 

S (p^Ts) 



Observing 

w{hK) — 1 



for h^O. 



Iih°' 

then, introducing the relationship of well-scaledness 

a{r,h) = ^^l (5.26) 
between the rescaling factors r and h, we finally get the limiting equation 



su(k, s) — 1 



\k\'^u{k,s), (5.27) 



corresponding to Eqs. (2.5) and (2.1) with 9 = 0, the symmetric space-time 
fractional diffusion equation. 
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Comments: Let us point out an advantage of splitting the passages 
r — > and /i — > 0. Whereas by the combined passage as in choice (C), 
if done in the weh-scaled way (5.26), the mystical concept of respeeding 
can be avoided, there arises the question of correct use of the continuity 
theorems of probability. There is one continuity theorem for the Laplace 
transform, one for the Fourier transform, see [8]. Possible doubts whether 
their simultaneous use is legitimate vanish by applying them in succession, as 
in our two splitting methods. For a more detailed discussion of mathematical 
aspects and the involved stochastic processes we refer to our recent paper 

m- 

6 Subordination in stochastic processes 

The common method for simulating particle trajectories consists in 
interpreting the concept of subordination (see [8j) as one of transforming 
a stochastic process y(t*) (we call it the parent process) where is not 
the physical but an operational time into a process X(t) in physical time 
t, by generating the operational time from the physical time t via a 
positively oriented stochastic process T^,{t) to arrive at the representation 
X{t) = Y{T^{t)). For simulation one then needs a routine for generating the 
process T,,(t). For simulating trajectories in space-time fractional diffusion 
(2.1) this requires simulation of the hitting time process which is inverse to 
the stable subordinator in Feller's parametrization [8j, the stable process of 
order (3 and skewness — /3. 

There are routines available for simulating stable variates, see e.g. [Ml [25]. 
But we do not know of easy routines for inverting a stable subordinator. 
Therefore, we have developed our method oi parametric subordination which, 
by starting from the operational time allows construction of trajectories by 
producing: FIRST an a-stable process x = Y{t^) with skewness 9 for for the 
position X in space, SECOND an extreme positive-oriented stable process 
t = T{t^) of order /? with skewness —(3 that we call the leading process. Then 
we get, in the (t, j;)-plane, the parametrized graph t = T{t:^),x = Y{t^) 
of a desired trajectory of the process X{t) corresponding to eq. (2.1) as 
X{t) = Y{T^{t)) where now t^ = T*(t) is the process 'inverse' to t = T(t,,). 
For the more general situation we refer to our recent paper [20]. This method 
is exact in the sense that it allows us to produce snapshots of a true particle 
path. Let us in this context also draw attention to the recent paper ^27j by 
Kleinhaus and Friedrich. 
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Let us sketch how this method directly arises from the CTRW approximation 
under appropriate power law assumptions on the waiting time and the 
jumps. Compare also with [39] and [17j. For mathematical details we refer 
to [20]. We start with the equations (3.4) and (3.5), and rescale waiting 
times and jumps again with positive factors r and h. In the Fourier-Laplace 
domain we assume, a bit more general than in our Master lemmata, the 
power-law conditions 

1 - 0(s) ~ a/ , A > as s ^ 0+ , (6.1) 

1-u}(k) ~/i[K|°i^^^SnK^ ^>Q^ as K^O. (6.2) 

We see that we have a combination of two Markov processes happening in 
discrete time, one giving a jump (with density w^x)) in space x at every 
instant n where n is the running index, the other one giving a positive jump 
in time t at every instant n. 

Using the effect of the rescaling on Eq. (3.5) and correspondingly decorating 
it with additional indices we get in the Fourier-Laplace domain 

%A^,s) = f; (^(^,))" (^)(/,^))« . (6.3) 

n=0 ^ 

Separately we treat the powers (J){ts)J and ({t;(/iK))", so avoiding 
the problematic simultaneous inversion of the diffusion limit from the 
Fourier-Laplace domain into the physical domain. Observing from (6.1) 
(^(piTs)^ ~ (l — A(rs)^)", we relate the running index n to the presumed 

operational time by n ~ t*/(Ar^), and for fixed s (as required by 
the continuity theorem of probability theory), by sending r ^ we get 

(j){Ts)j ~ (l — X s^y* ^ ^''''^ ^ — > exp (— t^, s^). Here the Laplace variable 
s corresponds to physical time t, and in Laplace inversion we must treat 
as a parameter. Hence, in physical time exp(— t,,s^) corresponds to 

5/3(t,i*) =C'^^^,0(C'^'^t), (6.4) 

with 5^(s) = exp(— s^). Here 5/3(t,t*) is the totally positively skewed stable 
density (with respect to the variable t) evolving in operational time 
according to the "space"- fractional equation 

J- gp{t, U) = tD^p gp{t, Q , ^^(t, 0) = 6{t) , (6.5) 
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where t plays the role of the spatial variable. Analogously, observing from 
(6.2) {w{hK))'' ~ (l- n{h\K\)''i^^^^^'^y, with the aim of obtaining a 

meaningful limit we set n ~ U/{iJ,h'^), and find, by sending h — > 0+, the 
relation 



{w{hK)r ~ (l ^exp(-t,|Kri^signK^ ^ 

the Fourier transform of a ^-skewed a-stable density fa,eix,t*) evolving in 
operational time i*. This density is the solution of the space-fractional 
equation 



_d_ 

3t* 



(6.6) 



The two relations between the running index n and the presumed operational 
time require the (asymptotic) scaling relation ~ that for 

purpose of computation we simplify to 



(6.7) 



Replacing U by = nXr^, using the asymptotic results obtained for 
the powers (j){Tsfj and {w{hK))^, furthermore noting {1 — (j)(Ts)) / s ~ 

■s^~^Ar'^, we finally obtain from (6.3) the Riemann sum (with increment 

At/3) 



n=0 



and hence the integral 



PhA'^^s) sl^' - I exp 
/o 



j exp — ^ 



^sign 



For the limiting process Ui3{x,t) this means 



Up{K,s) 



f 

Jo 



^ exp 



^sign 



(6.^ 



(6.9) 



(6.10) 



Observe that the RHS of this equation is just another way of writing the 
RHS of equation (2.6) which is the Fourier-Laplace solution of the space- 
time fractional diffusion equation (2.1). By inverting the transforms we get 
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after some manipulations (compare [31j ) in physical space-time the integral 
formula of subordination 

/•OD 

Ufi{x,t)=l fa,e{^^U)gp{t*,t)dt^ (6.11) 

JO 

with 

9p{U,t) = ygp {tt:''^) t:'/^-' . (6.12) 

Eq. (6.11) is basic for the conventional concept of subordination where there 
are also two processes involved. One is the unidirectional motion along the 
t^, axis representing the operational time. This motion happens in physical 
time t and the pdf for the operational time having value is (as density 
in t^, evolving in physical time t) given by (6.12). In fact, by substituting 
y = tt^ we find 

/ gp{t^,t)dt^= gp{t,t^)dt = l, Vt>0. (6.13) 
JO Jo 

The operational time t^ stands in analogy to the counting index n in Eqs. 
(3.5) and (6.3). The other process is the process described by Eq. (6.6), a 
spatial probability density for sojourn of the particle in point x evolving in 
operational time t*, namely 'U/3(x,t^=) = fa,e{x,t^). We get the solution to 
the Cauchy problem (2.1), namely the pdf u{x,t) = up{x^t) for sojourn in 
point X, evolving in physical time t, by averaging M^(x,t*) with the weight 
function gp{t^,t) over the interval < < oo according to (6.11). 

6.1 Trajectories for space-time fractional diffusion 

In the series representation (3.5) of the CTRW the running index n (the 
number of jumps having occurred up to physical time t) is a discrete 
operational time, proceeding in unit steps. To this index n corresponds the 
physical time t = tn, the sum of the first n waiting times, and in physical 
space the position x = Xn, the sum of the first n jumps, see Section 3. 

We have two discrete Markov processes (discrete in operational time n), 
namely a random walk in the space variable x, with jumps X„, and another 
random walk (only in positive direction) of the physical time t, making a 
forward jump T„ at every instant n. Rescaling space and physical time by 
factors h and r, observing the scaling relation (6.7), we introduce, by sending 
h ^ and r ^ 0, the continuous operational time 

t^^nXrl^ r^n^ih"' . (6.14) 
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Then in the diffusion hmit the spatial process becomes an a-stable process 
for the position x = x = x{t^:), whereas the unilateral time process becomes 
a unilateral (positively directed) /3-stable process for the physical time t = 
t = t{t^). A trajectory of a diffusing particle in physical coordinates can 
be produced by combining in the {t, x) plane the two Markovian random 
evolutions 

obeying the stochastic differential equations (compare with |27) ) 

dx = (i(Levy noise of order a and skewness 9) , , , 

dt = d[one sided positively oriented Levy noise of order (3) . 

This gives us in the {t, x) plane the t^, - parametrized particle trajectory that 
by elimination of t^^ we get as x = x(t). We suggest to call this procedure 
"construction of a particle trajectory by parametric subordination". Note 
that the process t = T(t,,) yielding the second random function in (6.16) 
has the properties of a subordinator in the sense of Definition 21.4 in |41j . 

Remark : It is instructive to see what happens for the limiting value (3 = 1. In 
this case the Laplace transform of gp{t, t^) = gi{t, t^) is exp(— t*s), implying 
gi{t,t^) = 6{t — t*), the delta density concentrated at i = i*. So, t = i*, 
operational time and physical time coincide. 



6.2 Numerical results for the symmetric case ^ = 

For numerical simulation of trajectories we proceed in three steps. 
First, let the operational time t^ assume N discrete equidistant values in a 
given interval [0, T], that is = nT/N, n = 0,1, . . . , N . As a working 
choice we take T = 1 and N = 10^. Then produce A'^ independent identically 
distributed (iid) random deviates, Yi,Y2, . . . , having a symmetric stable 
probability distribution of order a, see the book by Janicki [21] for a useful 
and efficient method to do that. Now, with the points 

n 

xo = 0, x„ = ^Xfc, n=l,...,A^, (6.17) 

k=l 

the couples {t^^n,Xn), plotted in the {t^,x) plane (operational time, physical 
space) can be considered as points of a true trajectory {x{t:^) : < t^ < 
T} of a symmetric Levy motion with order a corresponding to the integer 
values of operational time t^, = t^^n- In this identification of t^, with n 
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we use the fact that our stable laws for waiting times and jumps imply 
A = ^ = 1 in the asymptotics (6.1) and (6.2) and t = h = 1 as initial 
scaling factors in (6.3) and (6.7). In order to complete the trajectory we 
agree to connect every two successive points (t*,n)a^n) and x„+i) 
by a horizontal line from (t^^r^a^n) to a;„), and a vertical line from 

{t*,n+iTXn) to x„+i). Obviously, this is not the 'true' Levy motion 

from point (t*,n5a^n) to point (t*,n+i) a;„+i), but from the theory of CTRW 
we know this kind of discrete random process to converge in the appropriate 
sense to Levy motion. The points (t*,n) Xn) are points of a true Levy motion, 
as can be shown by observing the infinite divisibility and self-similarity of 
stable laws. 

As a second step, we produce N iid random deviates, Ti,T2, . . . , Tjv having 
a stable probability distribution with order f3 and skewness — /3 (extremal 
stable distributions). Then, consider the points 

n 

to = 0, tn = J2^k, n = l,...,N, (6.18) 

k=l 

and plot the couples (t*,n)in) in the (t*,t) (operational time, physical time) 
plane. By connecting points with horizontal and vertical lines we get 
snapsliots of trajectories {t(t*) : < < Nt =1} describing the evolution 
of the physical time t with increasing operational time t*. 

The final (third) step consists of plotting points (t(i*,n); 3;(i*^„)) in the (t, x) 
plane, namely the physical time-space plane, and connecting them as before. 
So one gets a discrete approximation of a particle trajectory of spatially 
symmetric {9 = 0) fractional diffusion with parameters a and p. 

Now as the successive values of t^, „ and generated by successively 

adding the relevant standardized stable random deviates, the obtained 
sets of points in the three coordinate planes: (t*,i), (t^,x), (t,x) can, 
in view of infinite divisibility and self-similarity of the stable probability 
distributions, be considered as snapshots of the corresponding true random 
processes occurring in continuous operational time and physical time t, 
correspondingly. Clearly, fine details between successive points are missing. 
They are hidden: 

- In the (t^,, x) plane in the horizontal lines from (t*,n) Xn) to (t*^„+i, Xn) and 
the vertical lines from (t^,^„+i, a;„) to (t*,n+i5 a;„+i). 

- In the (i*,i) plane in the horizontal lines from (t*,n)in) to (t*,n+i,in) and 
the vertical lines from (i*,n+i,in) to (t*^„+i, i„_|_i). 
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Figure 1: A trajectory for the parent process x = Y{t^,) with {a = 1.5}. 




Figure 2: A trajectory for the leading process t = T{t*). 
LEFT: {/3 = 0.90}, RIGHT: {/3 = 0.80}. 




Figure 3: A trajectory for the subordinated process x = X{t). 
LEFT: {a = 1.5, /3 = 0.90}, RIGHT: {a = 1.5, /3 = 0.80}. 
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- In the {t,x) plane in the horizontal lines from {tn,Xn) to {tn+i,Xn) and 
the vertical lines from {tn+i,Xn) to {tn+i,Xn+i)- 

The well-scaled passage to the diffusion limit here consists simply in 
regularly subdividing the {t*} intervals of length 1 into smaller and smaller 
subintervals (all of equal length r and adjusting the random increments 
of t and X according to the requirement of self-similarity, namely taking, 
respectively, the waiting times and spatial jumps as t^^^ multiplied by 
a standard extreme /3-stable deviate, r^/" multiplied by a standard (in 
our special case: symmetric) a-stable deviate, respectively, as required 
by the self-similarity properties of the stable probability distributions). 
Furthermore if we watch a trajectory in a large interval of operational time 
f*, the points (t*,^^^^) and (t*^„+i, a;„+i) will in the graphs appear very 
near to each other in operational time and aside from missing mutually 
cancelling jumps up and down (extremely near to each other) we have a 
good picture of the true processes. 

As interesting case studies let us present in Figs 1-3 the trajectories for 
{a = 1.5, P = 0.90,9 = 0} and {a = 1.5, /? = 0.80,6' = 0}, having in 
common the parent process x = Y(t^,). 

7 Conclusions 

Fractional diffusion processes as models of anomalous diffTision are gaining 
increasing popularity not only among science researchers but also among 
more or less pure mathematicians. For the latter they offer fascinating 
opportunities for applying pseudo-differential operators and other powerful 
analytic instruments, e.g. those of fractional calculus that in recent decades 
has made remarkable advances. 

In the field of anomalous diffusion we meet challenges for the experimental 
sciences, for mathematical modelling of real processes and their simulation, 
for investigation of the underlying evolution processes (the macroscopic 
aspect) and the fine-structure of their particle trajectories (the microscopic 
aspect), and for numerical analysis and computational treatment of less 
common problems. In our presentation we have discussed three topics of 
current interest that make visible the large arsenal of tools required. 
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Appendix: Proofs of the Lemmata 

We present here proofs for Master Lemma 1 and Master Lemma 2 of Section 
4. We had proved analogues of these Lemmata for the probability densities 
in [16]. Here we now give proofs for the probability distribution functions, 
thereby referring to Chapter 8 of the fundamental treatise by Bingham et 
al. [1] on Regular Variation through an appropriate change of notations. In 
view of this we first need to recall the notions of slowly varying functions 
and regularly varying functions. These concepts allow generalizations of the 
two Master Lemmata stated, without proofs, in [13] and |17j . 

Definitions: We call a (measurable) positive function o(y), defined in a 
right neighborhood of zero, slowly varying at zero if a{cy)/a{y) — > 1 with 
y ^ for every c > 0. We call a (measurable) positive function b{y), defined 
in a neighborhood of infinity, slowly varying at infinity if h{cy)/b{y) 1 with 
y ^ oo for every c > 0. An example of a slowly varying function at zero 
and infinity is: llogyp with 7 S R. Then regularly varying functions are 
power functions multiplied by slowly varying functions. 

Proof of Master Lemma 1: 

Note that because of symmetry we need only consider positive values of the 
variables x and n. 

In the easy case (a) q = 2 the well-known fact o"^ = — {^"(O) implies 
1 — w[K) ~ — /t as K ^ . 

In case (b) we refer to Theorem 8.1.10 in [3]. It says that if for a probability 
distribution function W(x) we set 

/+00 
cos(kx) dW{x) , 
-00 

and take any function L{x) slowly varying at infinity, then the relation 

T(x) ~ -Z^(a;) x^" for x ^ 00 , 
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is equivalent to the relation 

^ ' 2T{a) sin(a7r/2) ^ ' ' 

Taking now L{x) as the constant function L{x) = 2b /a and observing that 
because of our symmetry assumption on the jump distribution function 
W{x) we have W{—x) = 1 — W{x) and hence T{x) = 2[1 — 14^(x)] in all 
continuity points, we arrive at U{Hi) = w{k) and see that 

r+oo 

1 - W{x) = / dW{x') ~ ba-^x-"" for x ^ oo 

J X 

in view of aT{a) = T{a + 1) implies 

1 - id{K) ~ — — J— for K ^ 0+ . 

r(Q + 1) sm(Q;7r/2j 

We have completed the proof. 
Proof of Master Lemma 2. 

In the easy case (A) (3 = 1 the statement 1 — (j){s) = ps is a consequence of 
the well-known fact p = — 0'(O). 

In case (B) < /3 < 1 we invoke Corollary 8.1.7 of |3j. It says, among other 
things, that for a probability distribution function vanishing for t < 
the relation 

where L{t) is a slowly varying function at infinity, implies the relation 

1 - ^(s) ~ L{l/s) for s ^ 0+ . 
Now taking L{t) ee cr(l - /?)//? we get 

l-^(s)~A/ for s^0+, with A = c ^^-^~^^ = — ^^-^ttt , 

' ' 13 T{I3 + 1) sin(/?7r) ' 

where we have used the reflection formula for the gamma function. The 
proof is complete. 
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